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Abstract 



CLc We give a direct derivation of a proposal of Nekrasov-Shatashvili concerning the quantization 

^ ■ conditions of the Toda chain. The quantization conditions are formulated in terms of solutions 

j3 ' to a nonlinear integral equation similar to the ones coming from the thermodynamic Bethe 

ansatz. This is equivalent to extremizing a certain function called Yang's potential. It is shown 
that the Nekrasov-Shatashvili formulation of the quantization conditions follows from the 
^ ! solution theory of the Baxter equation, suggesting that this way of formulating the quantization 

O ■ conditions should indeed be applicable to large classes of quantized algebraically integrable 

ps) ■ models. 

O 

O . 

1 Introduction 

The A^-body quantum mechanical Hamiltonian 

e=l k=2 

is known as the quantum Toda chain. Above, and are quantum observables satisfying the 
canonical commutation relations [p^, Xk] = ihSi^k, 9 and k are coupling constant. When k = 1 
one deals with the so-called closed Toda chain and, when k = 0, with the open Toda chain. 

This model appears to be a prototype for an interesting class of integrable models called al- 
gebraically integrable models. It was introduced and solved, on the classical level, by Toda 



[Tlbl in 1967. Then, in 1977, Olshanetsky and Perelomov HOPII constructed the set of com- 
muting and independent integrals of motion for the closed chain, thus proving the so-called 
quantum integrability of the model. In 1980-81, Gutzwiller HGul was able to build explicitly 
the eigenfunctions and write down the quantization conditions for small numbers of particles 
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(A^ = 2, 3, 4). In particular he expressed the eigenfunctions of the closed chain with A^-sites as 
a linear combination of the eigenfunctions of the open chain with — 1 particles. 

A particularly succesful approach to the solution of the quantum Toda chain was initiated by 
Sklyanin. In 1985, Sklyanin applied the quantum inverse scattering method (QISM) to the study 
of the Toda chain. This led to the development of the so-called quantum separation of variables 
method. In this novel framework, he was able to obtain the Baxter equations for the model, 
from which Gutzwiller's quantization conditions can be obtained. The Baxter equation was 
later re-derived by Pasquier and Gaudin ffOP H with the help of an explicit construction of the 
so-called Q-operator, similar to the method developed by Baxter for the solution of the eight 
vertex model [|Bal . 

In '99 Kharchev and Lebedev nKL[| constructed the multiple integral representations for the 
eigenfunctions of the closed A^-particle Toda chain. Their construction can be seen as a gen- 
eralization of Gutzwiller's solution allowing one to express the eigenfunctions of the closed 
A^-particle Toda chain in terms of those of the open chain with — 1 particles for all A^. In '09, 
An HAnll completed the picture by proving rigorously that Gutzwiller's quantization conditions 
are necessary and sufficient for obtaining a state in the spectrum. 

However, the form of the quantization conditions obtained in the above-mentioned works ap- 
pears to be rather involved. Recently Nekrasov and Shatashvili proposed in HNS II that the quanti- 
zation conditions for the Toda chain can be reformulated in terms of the solutions to a nonlinear 
integral equation (NLIE) similar to the equations originating in the thermodynamic Bethe ansatz 
method. With the help of the solutions to the relevant nonlinear integral equation, Nekrasov and 
Shatashvili defined a function W whose critical points are in a one-to-one correspondence with 
the simultaneous eigenstates of the conserved quantitites. This formulation not only seems to be 
in some respects more efficient than the previous one, it also indicates an amazing universality 
of the form the quantization conditions may take in integrable models. 

The proposal of pNSII was based on rather indirect arguments coming from the study of su- 
persymmetric gauge theories. It seems desirable to derive the proposal more directly from the 
integrable structure of the model. Our main aim in this note is to give such a derivation. It 
is obtained from the solution theory of the Baxter equation. In other integrable models there 
are known connections between the Baxter equation and nonlinear integral equations that look 
similar to the one that appears here, see e.g. HBLZI |Zal |^. However, the precise form of the 
NLIE depends heavily on the analytic properties that the relevant solutions of the Baxter equa- 
tion must have in the different models. In the present case of a particle system we encounter 
an interesting new feature: the quantization conditions are not formulated as equations on the 
zeros of the solutions of the Baxter equation, but instead, they are equations on the poles of 
the so-called quantum Wronskian formed from two linearly independent solutions of the Bax- 
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ter equation. The positions 6 = {6i, . . . , 6^) of these poles are the variables that the Yang's 
potential W = W(5) depends on. 

It is worth stressing that the Baxter equation or generalizations thereof have a good chance to 
figure as a universal tool for the study of the spectrum of quantum integrable models. Our 
method of derivation strongly indicates that similar formulations of the quantization conditions 
should exist for large classes of quantized algebraically integrable models. 

This article is organized as follows. We first recall how the Separation of Variables method 
reduces the problem to find the eigenstates of the Toda chain to the problem to find a certain set 
of solutions to the Baxter functional equation specified by strong conditions on the analyticity 
and the asymptotics of its elements. Then, in Section [3] we explain how Gutzwiller's quantiza- 
tion conditions can be reformulated in terms of the solutions to a certain NLIE and in terms of 
the Yang's potential W(5). The derivation of this reformulation is sketched. The proofs of our 
claims are presented in the Appendices. In Appendix|Al we establish the relevant properties of 
Gutzwiller's basis of fundamental solutions to the T-Q equation. Then in Appendix|B]we prove 
the existence and uniqueness of solutions to the NLIE introduced in Section |3l In Appendix O 
we give rigorous proofs of the main results presented in Section [3l 

Acknowledgements. The authors gratefully acknowledge support from the EC by the Marie Curie 
Excellence Grant MEXT-CT-2006-042695. 

We are happy to dedicate this paper to T. Miwa on the occasion of his 60th birthday. 



2 Separation of variables approach to the Toda chain 
2.1 Integrability of the Toda chain 

The integrability of the Toda chain follows from the existence of Lax matrices 

L.(A)=( "T") [=^«'^'"]=^' (2.1) 

satisfying a Yang-Baxter equation with a rational, six-vertex type, R-matrix. Thus, the set of 
K-twisted monodromy matrices 

M(A)= J °, )^^(A)...Li(A) . (2.2) 

allows one to build the transfer matrix T (A) = tr [M (A)] which is the generating function of 
the set of commuting Hamiltonians associated with the Toda chain 

T(A) = A^ + ^(-l)'A^~^Hfc . (2.3) 

k=0 
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The first two Hamiltonians read 

N p2 f ^ 2 ^ 1 

Hi = EP'^ = P H2 = — - <^5^^ + + J2g2n^^.^r-^. (2.4) 

fc=l I e=l k=2 ) 

Any eigenvector of the transfer matrix defines a polynomial t (A) = Y\k=i ~ ''^fc)- ^ 
commuting Hamiltonians are self-adjoint, hence the set {r} is necessarily self conjugated: 

{Tk} = {Tk}- 

2.2 Separation of variables 

The Separation of Variables (SOV) method was developed for the Toda chain in llSkl IGP[ [KL[ 
I Anil . The main results of these works may be summarized as follows: 

The wave-functions \l/i ( ) of any eigenstate to the transfer matrix T(A) with 

eigenvalue t(A) can be represented by means of an integral transformation of the form 

^t(^) = / d/i(7)$i(7)Hp(7|x), (2.5) 

where integration is over vectors 7 = (71, . . . , '-)n-i) G M^^^ with respect to a measure d/i(7) 
first found in llSkll . Sp(7|x) is an integral kernel for which the explicit expression can be found 
in iKLjl, P is the eigenvalue of the center of mass momentum P in the state ^!^, and $f (7) is the 
wave-function in the so-called SOV-representation. The key feature of the SOV representation 
is that $^(7) takes a factorized form 

N-l 

$t(7) = n "itilk) . (2.6) 
fc=l 

The integral transformation (|2.5I) is constructed in such a way that the eigenvalue equation for 
the family of operators T(A) is equivalent to the fact that the function ^^(y) which represents 
the state $t via (12.51) and (12.61) satisfies the so-called Baxter equation, 

t(A)g,(A) = z^^?^\,(A + ih) + K\-ifg'^\{X - ih) . (2.7) 

The integral transformation (|2.5I) can be inverted to express $^(7) in terms of ^'^(x). In this 
way it becomes possible to find the necessary and sufficient conditions that g^d/) has to satisfy 
in order to represent an eigenstate of T(A) via (12.51) and (12.61) . The conditions are HKLl lAnH 

(i) t(A) is a polynomial of the form t{\) = Y[k=i('^ ~^ '^k), with {t^} = {fk}- 

(ii) g(A) is entire and has asymptotic behavior |g (A) I = 0(e~^l'^(^)l |_x|f (2|SWI-'^) 

Above, the O symbol is uniform in the strip {z : {z)\ < h/2}. This reduces the problem to 
construct all the eigenfunctions and finding the complete spectrum of the Toda chain to finding 
the set S of all solutions {t (A) , (A)) to the Baxter equation (|2.7I) that satisfy the conditions 
(i) and (ii) above. 
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3 Quantization conditions 



It turns out that the Baxter equation (12.71) admits solutions within the class § described above 
only for a discrete set of choices for the polyonomial t(\). This is what expresses the quanti- 
zation of the spectrum of T(A) within the SOV-framework. Our first aim in this section will be 
to outline how to reformulate the resulting conditions on t{\) more concretely, following the 
approaches initiated by Gutzwiller and Pasquier-Gaudin. 

It gives useful insight to divide the problem to construct and classify the solutions to the Baxter 
equation (12.71) which satisfy (i) and (ii) into in two steps. In the first step, one weakens the 
analytic requirements (ii) on g(A) slightly by allowing g(A) to have a certain number of poles. 
In this case, it will be possible to find two linearly independent solutions q^{X) to (|2.7I) for 
arbitrary t(A) satisfying (i). In the second step, one constructs the solution g(A) satisfying (i) 
and (ii) in the form 

g(A) = P+g+(A) + P_g-(A), (3.1) 

where P± are constants. The requirement that g(A) is entire means that the poles of q^{X) must 
cancel each other in (13.11) which is only possible if t(A) is fine-tuned in a suitable way. This is 
the origin of the quantization of the spectrum of T(A). 

3.1 Gutzwiller 's formulation of the quantization conditions 

It turns out that there is a canonical minmal choice for the set of poles of (A) that one needs 
to allow. One needs to allow poles 5i, . . . ,5n whose positions are determined by the choice 
of t(A). More precisely, out of t(A) one constructs the so-called Hill determinant 



n (A) = det 



P" 



t{X- ih) 




t(A) 



1 



t (A - ih) 
1 




1 



t(A) 



0... 



(3.2) 



where p := ng . It can be shown that ^(A) admits the representation 



TT 



N sinh - (A - (5a 

^(A)=n- 



TT 



a=i sinh -^{\- To) 



(3.3) 



where one chooses (5^)1 < fi/2. This defines 5i, . . . ,5n in terms of t{\). 
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Let us then, instead of S consider the class S' of solutions to (12.71) which satisfy the conditions 
(i) and (ii)', 



(ii)' g(A) is meromorphic with set of poles contained in . . . , 5^} and 



it has an asymptotic behavior \q (A)| = 0(e '2ft |A| 



-{2\^(X)\-h) 



)• 



The Baxter equation (12.71 ) has two linearly independent solutions qf{\) within S' for arbitrary 
t{\). One possible construction of the solutions qf{X) goes back to Gutzwiller's work on the 
Toda chain. They may be defined as follows: 



nL{e-^sinhf(A-5.)} 



(3.4) 



where 



n h-'T.r{i-t{X-rk)/h) 

k=l 



giNX^_ (A) e 



h 



n h'r.T{i + t{X-n)/h) 

k=l 



(3.5) 



with K± (A) being half-infinite determinants: 



K+ (A) = det 



P" 



t{X + 2ih) 




t-^{X + ih) 

1 t-^ (A + 2ih) 



(3.6) 



andK^ (A) = (A) (recall that {r^} = {r^}). For the reader's convenience we have included 
a self-contained proof that qf G S' in Appendix A. It's worth noting that Qf are linearly 
independent entire functions whose Wronskian can be evaluated explicitly, cf. Lemma [TJ 

W [Qt, g,1 (A) = K-^'g-'^h-'''^^' • n { ^ ^iX-n).n (A) | , (3.7) 

a=l ^ 



It can be shown that the most general solution q{X) G S' to the Baxter equation (12.71) may be 
represented in the form (13.11 ). The additional requirement that g(A) should be entire implies 



Qt{6a)-CQt{6a) = 0, for a = l,...,iV and some C e C, |C| = 1 



(3.8) 



to be supplemented by the condition that 6k = P HGui IGPH . This formulation of the 

quantization conditions looks fairly involved. It may be considered as a highly transcendental 
system of equations on the parameters determining t(A), in which both Q^{X) and 61, . . . ,6n 
have to be constructed from t(A) by means of (13.51) and (13.21) . (13.31) respectively. 
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3.2 Reformulation in terms of solutions to a nonlinear integral equation 



One may note that the set of parameters in 6 = (5i , . . . , ^a?) is just as big as the set of parameters 
in r = (ti, . . . , r^v) characterizing the polynomial t(A) appearing on the left hand side of the 
Baxter equation. The form of the quantization conditions (13.81) suggests that it may be useful 
to formulate these conditions directly in terms of the parameters S = (Si, . . . , S^) with t(\) = 
t{\\T (S)) being determined in terms of S by inverting the relation S = 5{t). A more convenient 
representation of the quantization conditions (13.81) would then be obtained if one was able to 
construct the solutions Qf{\) more directly as functions of the parameters S = (Si, ... , 5^)- 
In the following, for a given polynomial ?9(A) = nf=i('^ ~ ^fc) "^i^h complex conjugated roots, 
we will construct functions Qf{\). These will be shown to yield solutions the Baxter equation 
(12.71) via (|3.4I) . with ^^(A) being a polynomial whose coefficients depend on the parameters 5. 

The functions Qf{X) will be build out of the solutions ^^(A) to the following NLIE, 

log UX) = I m - ,) 1„ (l + . (3.9) 

where 

It will be shown in Appendix |B] that the solutions 15(A) to (13.91) are unique, and that they exist 
for all tuples 5 = {5i, . . . , 5^) of zeros of Hill determinants 'H(A) constructed from polynomials 
t(A) whose zeroes satisfy |53 (r^)! < h/2. The function 15(A) is meromorphic, with its poles 
accumulating in the direction |arg (A)| = n/2 and such that — )• 1 if A — )• 00 for A uniformly 
away from its set of poles. The properties of I5 allow one to define two auxiliary functions: 



ln.,(A) = - '---(1 + 

' J 2i7cX-fi + ih/2 V 



2i7r A - /i + ih/2 \ ^9 (/i - ih/2) ^9 (/i + ih/2) 
Invi (A -th) = / — — - 1 + 



2i7r A - /X - ih/2 \ d - ih/2) + ih/2) 



(3.11) 



Out of f| (A) and (A — ih), we may then construct 

n r (1 - 2 (A - 4) /^) n r (1 + ^ (A - 5fe) /h) 

k=l k=l 

It is shown in Appendix O that the functions Q"^ are entire (c/. LemmaH)). It is also shown there 
that the functions qf{X) defined from Q^{X) by relations like (13.41) are solutions to the Baxter 
equation (12.71) which have the right asymptotic behavior to be contained in S'. 
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One may therefore take the construction of Qf{X) from the solutions of the nonlinear inte- 
gral equation (13.91) as a replacement for the construction based on Gutzwillers solutions. The 
quantization conditions (|3.8I) may now be rewritten in terms of ^^(A) in the form 

2vrn. = ^ln/l-41np + zlnC-zX.l^ r(l-z(4-^,)M) ^^"'^^ 



as is fully demonstrated in Appendix O This form of the quantization condition may be more 
convenient for many applications than the ones previously obtained, equations (13.81) . 



3.3 Solutions to the Baxter equation from the solutions to a NLIE 

In order to understand how the connection between nonlinear integral equations and the Baxter 
equation comes about, the key observation is that the two functions qf defined above constitute 
a system of two linearly independent solutions of the Baxter equation (12.71) . This fact can be 
deduced from the so-called quantum Wronskian equation satisfied by Qf: 

/J- _2iA \NN 

Qt (A) QJ (A + th) - Qs (A) Qj (A + th) = K-'^ ( -^—j- ] H sinh ^ (A - 6^) . (3.14) 

The quantum Wronskian equation allows one to show that Qf{\) satisfy a Baxter-type equation 

tsWQfW = ^-''g'''Qf{X + ih) + K\i)''g'''^Q^{\ - ih) , (3.15) 
where the polynomial (A) is defined by 

, _ Qs (A - ifi) Qj (A + th) - Qj (A + ih) Qj (A - ih) 

' ^ ^ ~ ^ ' ^ Qt (A) QJ (A + ^h) - Qj (A + ^h) Qj (A) " ^'-''^ 

On the one hand, the Wronskian relation (13.141) allows one to show that ^^(A) and Qf{X) are 
related by the Baxter equation (13.151) . On the other hand, it also ensures that the residues of 
the possible poles of ts (13.161) vanish. Then, the polynomiality of ts is a consequence of the 
asymptotic behavior of the functions Qj. The details of the arguments are found in Appendix 

o 

In order to see how the quantum Wronskian relation is connected to the NLIE (13.91) . let us, 
starting from a solution Ys{X) to (13.91) . introduce two functions ft(A) and f4,(A) via (13.111) . On 
the one hand, noting that the kernel K(X) defined in (|3.10l) can be written as 

K{X) ^ 



2ni \X — ih X + ih 
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it is easy to see that (13.91) implies 

In Ys (A) = In [y^ (A + ih/2)) + In {v^ (A - 
On the other hand, note that 



In 



t^t ( A - + iO 



+ In 



-iO 



n (a - 4 

d/i 1 
2iTT A — 



(3.17) 



(3.18) 



P'n (p) 



?9 (/i - in/2) ^{p + ih/2) 



1 + 



(A) 



t9 (A - iV2) ?9 (A + 

Thus, using that v^j^ are meromorphic on C, we are able to continue the obtained relation 
everywhere on C , leading to the functional relation 



(A) (A) = 1 + 



P 



-v^ (A + ih) vi (A — ih) 



(3.19) 



Rewriting this in terms of Q"^ by means of (|3.12l) yields the quantum Wronskian equation (13.141) . 

At the moment we don't have a direct proof that a solution to (13.91 ) exists for all choices of ^9(A). 
We are able, however, to prove that all functions ^^(A) that can be constructed from Gutzwillers 
solutions are in fact solutions to (|3.9I) . This implies that all functions 15(A) needed for the 
formulation of the quantization conditions (13.131 ) can be obtained in this way. 

Let us finally note that there is a more direct way (c/. Proposition [5]) to reconstruct the Newton 
polynomials in the zeroes {rfc} of ts from the solution Ys to (|2.7I) : 

N N „ , 

P'n (r) 



(3.20) 



X In 1 + 



\^{T-ih/2)[ 



As one may reconstruct the eigenvalues of the conserved quantities from the Yl!^=i 
this essentially amounts to a reconstruction of the hk. 



3.4 Definition of Yang's potential 

It is interesting to notice that the quantization conditions (13.131) characterize the extrema of 
a certain function W (5) called Yang's potential in UNSL This Yang's potential is defined as 

W (S) = Vy'°'^* (6) + VyP"''* (6), where 

W^^''{5) = iY,^^^[ ] -lnC$^4+ 5ZG7(4-5,)-2z7r^nfc, (3.21) 

k=l VP/ ■ f,^^ 
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where w' (A) = In F (1 + i\/h), rik are some integers paremeterizing the eigenstate, and 



W^'^"* (5) = (3.22) 
lny^(/i) L ^ P^'Ys (/i) \ ^ ^.^ ^ -p'-^r^ (/i) A 1 d/i 



|^?(/i + zn/2)|V V|^9(Ai + ^V2)|V J 2ivr 



where 



?ln(l-t) 

Li2 (2) = / ^dt . (3.23) 

It seems worth emphasizing that, in our treatment, the reformulation of the quantization condi- 
tions in terms of Yang's potential was obtained from the solution theory of the Baxter equation 
(12.71) . This may be seen as a hint towards a more direct understanding of the claim in [INSll that 
the quantization conditions for large classes of quantized algebraically integrable models can 
be formulated in this way. The claim should follow quite generally from the solution theory of 
the Baxter equation. 



A Properties of Gutz wilier 's solutions 
A.l Analytic properties of Gutzwiller's solution 

The explicit construction of a set of two linearly independent entire solutions Qf of (|2.7I) with 
arbitrary monic polynomial t (A) goes back to Gutzwiller HGull . Prior to writing down these two 
solutions, recall the definition of the Wronskian of two solutions qi and 52 defined by 

W [gi, gs] (A) = gi (A) (A + ih) - gs (A) gi (A + ih) . (A.l) 

It is straightforward to see, using (12.71) . that W [gi, g2] is ih quasi-periodic: 

W [gi, ga] (A + th) = {-i f ^^W [q,, q^] (A) . (A.2) 

Proposition 1. Let t (A) = Y[k=i (-^ ^ ''"fe) monic polynomial of degree N with roots ap- 
pearing in complex-conjugate pairs {rfc} = {fk}. Then, the two functions below are entire 
solutions to the Baxter equation d3.i5l) . 

n h-'iT{i - z(A - Tk)/h) n ^'^r(i + z(A - n)/h) 

k=l k=l 
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Here K± (A) correspond to the unique meromorphic solutions to difference equations 

go to 1, when A — oo uniformly away from their set of poles. The solutions to these recur- 
rence relations are given explicitly by the determinant formula (13.61) and its complex conjugate. 
These two linearly independent solutions are entire and posses the asymptotic behavior 

|g±(A)| =e-'^«WO(e+^l«WI|A|^(^'^(')-")) for 3? (A) ^ ±00 , (A.6) 

and the O is uniform in every bounded strip ofC 

Proof. The only non-trivial part concerns the asymptotic behavior of K±. It follows as a corol- 
lary of Lemma [3l □ 

Lemma 1. The functions Qf are linearly independent and their Wronskian is expressed in 
terms of the Hill determinant di.2l) by (I3.7|) . The latter is closely related to and K_ 

Its zeroes form complex conjugated pairs {6k} = {Sk}, belong to the fundamental strip 
{z : {z)\ < h/2} and fulfill Xl^^i = XlJLi V 

Proof. The fact that Qf are linearly independent is a consequence of the fact that their Wron- 
skian does not vanish identically. The explicit expression for this Wronskian follows after some 
algebra. 

As we have assumed that the set {rfc} is self-conjugated, it follows from the determinant repre- 
sentation for "H that T-L (A) = "H (A), ie, the set {5k} is self-conjugated. In its turn, this implies 
a particular relation between the set of 5's and r's. Namely, computing the 3? (A) — )■ +00 
asymptotics of "H yields 

N N 

Tp = 6p + inh , for some n G N . (A.8) 
p=i p=i 

However, as X] ''"fc ^ K and ^ 5^ G M, the only possibility is n = . □ 
We are now in position to prove the 

Lemma 2. Let q be any meromorphic solution to (12.71) . Then, there exists two meromorphic 
ih-periodic functions P± (A) such that 

q (A) = P+ (A) Qt (A) + P (A) g,- (A) . (A.9) 
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Proof. Let q be any meromorphic solution to Baxter's T-Q equation (12.71) . Then consider 

»'[«,e.1(A) ' w[Qt,Q;][\) ^'^ ' 

The ratio of two Wronskian being ih periodic, one gets that, by construction 

W [g, Qt] (A) = W [g, Q^] (A) = . (A.ll) 
This leads to the system of equations for g (A): 



g+(A) g+(A + z;i) \ / -q{\ + ih) 
g-(A) g,-(A + z;i) i I g(A) 



(A. 12) 



Given any fixed A, there exist non-trivial solutions to (IA.12I) if only if the determinant of the 
matrix defining the system vanishes, ie W [g^, Q^^^ (A) = 0. However, it follows from (13.71 ) 
that W [g^, g^] (A) is an entire function that is non-identically zero. Therefore, it can only 
vanish at isolated points. Hence, we get that g (A) ^ only at an at most countable set. As g (A) 
is meromorphic on C, g = 0. □ 

We now provide a rough characterization of the set of zeroes of Qf. As the V function has no 
zeroes on C, the only zeroes of Qf are those of K± (A). 

Proposition 2. Assume that (r^)! < h/2 and that the set {r^} is invariant under com- 
plex conjugation. Then, the set of zeroes of (A — ih/2) belongs to the half-plane 

[zeC : '^{z) < -h/2} and 

|/f,(A-«/2)r^^_ for A.R. (A.13) 



H{X-th/2) 

A similar statement holds for K_ {X -\- ih/2), namely the set of zeroes ofK_ {X + ih/2) lies in 
the half-plane {z & C : Q (z) > +h/2} and (A + ih/2) does not vanish on M. 

Proof. It follows from the determinant representations that (X — ih/2) has poles at — 
i {2n + 1) h/2, n G N, in particular they all belong to the half-plane {z E C : Q (z) < —h/2}. 
Also, since the zeroes and poles of the Hill determinant are self-conjugated, 

N cosh ^ (A -4) 

'H{X-ih/2) = Yl w >0 ' VAgM. (A.14) 

fc=i cosh - (A - Tfc) 

As the set {rk} is self-conjugate, it is easy to see that 



K+{X) = K^{X) and t {X - ih/2) t {X + ih/2) = \t {X - ih/2)\\ (A.15) 



12 



This allows us to rewrite (IA.7I ) in the form 



\K^{X-ih/2)[ 
H{\-ih/2) 



P" 



n{x-ih/2) 



K+ (A + ih/2) 



t (A - ih/2) 



(A.16) 



It follows from (IA.14I ) and (IA.16I) that there exists a c > such that |K+ (A - i/i/2)| > c for 
A G R. Hence, (A — has no zeroes on M. As Kj^ (A — ih/2) has manifestly no poles 
in the half-plane {z : ^ [z) > —h/2}, one has that 



/(P) 



(A.17) 



Here, we remind that p is the deformation parameter appearing in (|3.6I) . We also specify that 
the function / (p) is well defined as \K^ 



+ 1 |K-tS/2 



> c and the ratio K'j^/ decays at least as 
j^-{2Af+i) ^j- infinity, uniformly in p, cf lemma [3l Thus, applying the dominated convergence 
theorem we obtain that / (p) is continuous in p. As it is integer valued, it is constant. The value 
of this constant is fixed from / (0) = (as then Kj^ = 1). This shows that (A) has all of its 
zeros lying below the line M — □ 



A.2 Bounds for 

Lemma 3. Let {rk} = {rk} and (r^)! < h/2, then K± ^ for A — oo uniformly away 
from its set of poles and 

g (A) = O (A-(^^+^)) , (A.18) 
where the O is uniform as long as p belongs to some fixed compact subset ofC 

Proof. As the T-Q equations can be solved explicitly when = 1 in terms of Bessel functions, 
it is enough to consider the case N > 2. We focus on as behavior of K_ follows by complex 
conjugation. 

admits the discrete Fredholm series representation: 

if + ( A) = 1 + J] ^ J2 det„ [M,,,, ( A)] , (A. 19) 



ni 

n>l hi,...,h„eN 



where we have 



A = . + — r — . (A.20) 

c (A — lah) r (A — lah) 
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Then, by Haddamard's inequality 



J2 det„ [Mh^h, (A)] 

hi,...,h„eN 



n I n 

s E n Ei^f' 

/ii,...,h„GNa=l L 6=1 
n ( n 



haht 



a=l \^a,b=l 

There we have set 

+00 / (-> \ A'' +00 



(A) = J]|t(A-iA;;i)p^ < ( - j 5^A;-^ and u (A) = O (^A^^) . (A.22) 

A;=l ^ ^ k=l 



Hence, we get that 

\K+ (A) - 1| < e"(^Ki+l^l') - 1 , (A.23) 

and thus — )• 1 for A — t- cxo uniformly away from its set of poles. It remains to provide the 
stronger estimates for its decay at infinity. 

Termwise differentiation of the series leads to 

K w I < E S (^) (1 + \p\T") ^ w (1 + IpH < s (A) (1 + /) e-w(^+i^i") , 

n>l 

where u (A) = J2n>i 1^'/^^ (A — "inh) | = O (|Ap^) . One then takes the derivative of the Hill's 
determinant relation (IA.7I) at X — ih/2. This leads to 

The uniform estimates in p that we have established combined with the fact that (A) = 
O (A~°°) uniformly in p lead to the desired form of the estimates, with a O that is uniform as 
long as p belongs to some compact subset of C. □ 



B Existence and uniqueness of solutions to the nonlinear in- 
tegral equation 

In this Appendix we prove the existence and uniqueness of solutions to the TBA non-linear 
integral equation (|3.15l) . Let 1!) (A) = 11^=1 (A — ^fc) have its zeroes given by the N zeroes of 
the Hill determinant built out of i(A). Then set 

n{x-tn/2) t{\-ih/2)t{x + ih/2) ' ^ ^ ^ 
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It is a straightforward consequence of Lemma [3] that Yt is a meromorphic function whose poles 
accumulate in the direction |arg (A) | = n/2 and such that It — )• 1 for A — )■ oo uniformly aways 
from the set of its poles. 

Proposition 3. The function In Yt defined in ^^.71) is continuous, positive and bounded on M. It 
is the unique solution in this class to the non-linear integral equation ^3.9P . 



Proof. We first prove the uniqueness of solutions. Let || ■ ||^ stand for the sup norm on bounded 
and continuous functions on M. We set J" = {/ G (M) : / > and ||/||^ < +oo}. 
Then we define the operator L on by 



L [/] (A) 



yd,;.-(A-,).n(l. j. (B.2) 



The mapping L stabilizes J^. Indeed, 



d/.ir(A-/.)^^^^— ^<ln(l + pM/II^J-i) , (B.3) 



where J = infxeR W - ih/2)f > 0, due to |^ {5k)\ < h/2. 
Any solution to the NLIE appears as a fixed point of L in J^. We shall now prove that L can 



have at most one fixed point. This settles the question of uniqueness of solutions to (13.91) . This 
part goes as in HFKSH . Let f,g E then 



\L[f]-L[g]\{\) 



^ p^e9{r)+t{f-g){T) 



dt jdrK (A - r) . , , (/ - g) (r) 







\d (r + in/2)y + p^e9{r)+t{f-g){r) 

max( IIqII^,,!! f II 
J + p/l,emax(||5lU,|| 



r.max(||5||^, 11/11^) 

< ^ . „.„ , \\f-9\\oo<\\f-9\\^- (B.4) 



Hence, L admits a unique fixed point. 

A direct proof of existence of the solutions to (13.91 ) is possible if p'^/J < 1. In this case it is 
easily seen that L [/] (A) is a bounded mapping in the sense that it stabilizes all balls in F of 
radius R> — In (l — p^' / j) . In such a case, (IB. 41) implies that L [f] is a contractive map on a 
Banach space. It thus admits a unique fixed point. 

However, it is always possible to construct a solution to (13.91) in terms of the the half-infinite 
determinants K±. Recall that {rk} and hence {6k} are invariant under complex conjugation. 
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Let 



k=l 

N 



i (A - Tk) /h) 

r(2 (A - Sk) /h) 



It follows from Proposition [2] that v^/i (A — ih/2) are holomorphic and non-vanishing in ]HI+/_. 
Moreover, as ^ 5^ = ^ r^, we get that f-f/j, (A — ih/2) = 1 + (A^^) in their respective 
domains of holomorphy. Also, due to the Hill determinant identity (IA.7I) 

We agree upon choosing a determination of v^/^ (A — ih/2) such that 

Inv-f/l (A - ih/2) — ^ ^ In [v-^v^] (A - ih/2) = Inv-, (A - ih/2) + Inv^ (A - ih/2) . 

A— S-+00 

Thus, for A G M, by computing the residues in the upper or lower half-plane and using the decay 
properties of the integrand at infinity, one sees that and are recovered from Yt{\) via the 
integral representations (13.1 II) . Hence, with the same choice of branches of logarithm as before 
(the one that goes to when 9? (A) goes to +oo) we get, on the one hand, that 

In K (A + th/2) (A - 3zh/2)] = [ drK (A - r) In (l + ^^XlizL_\ . (b.8) 

J V \i9{T-th 2)\ J 

R 

On the other hand, it is straightforward to check that (A + ih/2) (A — 'iih/2) = Y (X) . 
This proves the existence of the relevant set of solutions to (13.91) . □ 



C Baxter Equation and quantization conditions from TBA 

We first prove basic properties of the functions Qf. Then we derive the T-Q equation generated 
by Qf and finally obtain the quantization conditons. 

C.l Analytic properties ofQf 

Lemma 4. The functions Qf defined in (13.121) are entire and have the asymptotic behavior 

\Qf\=e-'^ ■ 0(e^t«WI|A|^(±^^(^)-'^)) 5?(A)^±oo (C.l) 
where the O symbol is uniform in {z : 1^3 (2;) | < h/2}. 
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Proof. It is readily seen from the asymptotic behavior in the strip {z : \'^{z)\ < K] of the 
solution Ys to (13.91) . that (A) -> 1 and v^^ (A — ih) — )• 1 when 3? (A) — )■ ±oo in the strip 
{z : {z)\ < h/2}. Then a straightforward computation leads to (|C.1I) . We assume that all 
the (5's are distinct and, if necessary, take the limit of coinciding 5's at the end of the calculation. 

It remains to prove that Qf are entire. For this, we show that the products of F-functions cancel 
the poles of v^/i. We need to construct a meromorphic continuation to C of v^/i starting from 
the strip B^^K with = {z : \^{z)\<nh}. It follows from the very form of the NLIE (Il9l) 
that the solution Yg is holomorphic in B^^\ Let us introduce the notation 

The only singularities of Vg (/i) /Vs (jf) in B^-^^ correspond to the zeroes of Vs and to its poles. 
The latter are located at = 5^ ± ih/2, k = 1, . . . , N. Hence, as Vg (fi) /Vs (yu) is decaying 
sufficiently fast at infinity, one gets 

Yg , ^ sr^ n,, sr^ n 



Ys ^ X — z — ih ^ \ — z + ih 

Vs{z)=0 Vs{z)=0 

n N ^ n N 

"^1? A-4-^(2p+l)/l/2"^^A-4 + ^(2p+l)/l/2 



1 



+ 



d/i V^ (/i) /Vs ifi) f d/i V^ (/i) /Vs (/i) 



2in \ — n — ih J 2iTi \ — fi + ih 



Above, we have denoted by riz the multiplicity of a zero ^ of and B^\ resp. B^^\ stands for 
n e+, resp. n H_. It thus follows that Vs has simple poles at 4 + if (2n + 1), n G Z. 
Exactly the same reasoning as before shows that v-^ (A — ih) has its only simple poles at 5k—inh, 
k = 1, . . . , N and n eW. Similarly, (A) has its only simple poles at 6k + inh, k = 1, . . . , N 
and n e'N. Therefore, these poles are canceled out by the zeroes of the F-functions and Qf are 
both entire. □ 



C.2 Baxter equation 

Proposition 4. The polynomial ts (A) given in (13.161) is a monic real valued polynomial of 
degree N. It is such that the functions Qf solve the Baxter equation 

ts{\)Qf{\) = ^''g''"Qf{\ + ^h) + ^\-^fg''"Qf{\-^h), (C.4) 

Finally, given any monic real valued polynomial of degree N with roots in the strip 
{z : |53 {z)\ < h/2}, it is always possible to find a complex-conjugation invariant set of pa- 
rameters {6k} such that ts equals to this polynomial. 
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Proof. We first prove that Q'l satisfies (IC.4I) with ts being given by (13.161) . 

n I - sinh ^ (A - 5p) I ■ U (A) Qj (A) (C.5) 



AT 



p=l 



_22rA \ ^ AT 



(/s:/^)" g+ (A - ih) { I j n^inh ^ (A - 5,) + g+ (A + ^/i) g,- (A) 



iV 



^ ^-e . •■J||^sinh^(A-(5p)|{(7^'"'z^g+(A + z;i) + /€''(7^'^(-0^g+^ 

In the intermediate steps, we have used the definition of and the q-Wronskian equation (13.141 ). 

Now we show that t^, as defined in (13.161) . is indeed a monic real valued polynomial of degree 
A^. We first assume that the are pairwise distinct. The case when several bk% coincide 
follows by taking the limit in the final formulae. As are both entire, we get that the only 
potential poles of (A) are located at A = + ^nfi, n ^TL. The set of zeroes of v^j^ differs 
necessarily from its set of poles (as follows readily from (IC.3I) and similar representations for 
v^l^. Hence, {dk + inh) ^ 0, for k = 1, . . . , N and n E Z. Therefore, it follows from the 
Wronskian relation (13.141) , that 

TwJT = I • fo^^ /c G 1 1 ; iV] and neN . (C.6) 

This implies that the possible poles of the expression in (13.161) get canceled. It follows that ts is 
entire. It remains to control its asymptotic behavior. We may express ts in terms of v^/^, 

ts (A) = v^{\- th) (A) n (A - 5a) - ■ 

^^=1 n (A - 5a) {h^ + (A - 6af) 

a=l 

Due to the asymptotic behvavior of Ys dXoo, one can deform the integration contour in the def- 
inition v^/i so as to obtain its asymptotic behavior in the whole plane A — )■ oo, for A uniformly 
away from the set of poles of v^^. 

As we have that f^/^ — t- 1 when A — )■ oo, we get that ts — \^ when A — )■ oo. Hence, ts is a 
monic polynomial of degree A^. It is real valued for A G M as for such A's, Vq (A) G M, what 
implies that (A — ih) = vi^ (A), ie ts (A) = ts (A) . 

The fact that, in this way, one is able to generate any monic polynomial with roots in the 
strip {z : |53 (2;)| < h/2] follows form the uniqueness of solutions to the TBA-NLIE and the 
construction of the function v^^^ in terms of determinants, as given in (IB.5I) - (IB.6I) . □ 
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C.3 The quantization conditions 

We now prove that the quantization conditions for the model (conditions on the zeroes of the 
polynomial ts (X) for (12.71) to have entire solutions with a prescribed decay as given in point (ii)) 
can be written down in a TBA-like form. Moreover, as opposed to the Gutzwiller form of the 
quantization conditions, the ones that will follow only involve one set of parameters. Namely, 
the zeroes {5k} of the Hill determinant associated with ts (A) given in (13.161) . We show that 
under certain reasonable assumptions, it is possible to reconstruct the Newton polynomials in 
the zeroes of ts and hence the spectrum of the model. This proves the Nekrasov-Shatashvili 
conjecture UNSH . We first reconstruct the zeroes of ts. 

Proposition 5. Let ts (A) = H^i (-^ ~ '^a) ^ polynomial whose zeroes lie in the strip 
{z E C : 1^3 (2;) I < h/2}. If{Sk} is the associated set of zeroes of the Hill determinant and Ys 
the unique solution to the NLIE (13.91 ). then the Newton polynomials 8k = J2p=i '^p zeroes 
ofts are reconstructed by means of formula ( 13.201) above. The convergence of these integrals is 
part of the conclusion. 



Proof. Due to the uniqueness of solutions to the NLIE (13.91 ). one has that the solution Ys can 
be expressed, as in (IB. II) in terms of K±, T-L. The latter determinants are parameterized by the 
zeroes {rk} of ts, and the parameters {5k} appearing in the NLIE (13.91 ) coincide with the set of 
zeroes of the Hill determinant. By invoking the continuity of the logarithm on R and its decay 
at infinity, we get 



d/i 
2i7r 



-ih/2 



(C.8) 



R+ih/2 



In the intermediate steps, we have used the quick decay at infinity of the integrand 

g (A) = O (A-^-^) and ^ (A) = O (A^) . 



(C.9) 



This allows us to split the integral in three and compute the parts involving K+, resp. K^, by 
the residues in the upper/lower half plane (thus giving 0). Hence, the only part that gives a non- 
trivial contribution is the contour integral involving H' /H. The only poles that contribute to the 
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result are located at the zeroes 5^ of the Hill determinant (they have residue +1) and at the poles 
Tfc of the Hill determinant (they have residue -1) that are located in the strip {z) \ < fi/2. □ 

This result offers a direct way to recover the spectum of the model from a solution to the TBA 
equation (13.91) . It remains to derive the set of quantization conditions on the parameters 5k- 

Theorem 1. There exists a unique entire solution q to the T-Q equation (|2.7I) whose asymptotic 
behavior is as stated in (ii) if and only if the parameters {6 k} appearing in the TBA NLIE (|3.9I) 
satisfy to the quantization conditions given in (13.131) . 

Remark 1. The solvability of the quantization conditions, the occurrence of complex solutions 
( 3= {5k) 7^ 0, G {-2 : |Q= (-2)1 < h/2}, the uniqueness of solutions for a given choice of 
integers G Z are all open questions. 

Proof. According to lemma l2l any meromorphic solution q to the T-Q equation takes the form 

" - W[Qt,Q^] (A) * W[«,Q7](A) * ■ '"-'"^ 

Recall that the Wronskian W \_Q'l , Q^] (A) is given by (|3.14l) . It is possible to compute the 
Wronskians W [g, Q^] (A) by using the asymptotic behavior of q and Qf. Due to their ih 
quasi-periodicity, these Wronskians take the form W [g, Qf]^ (A) = e~^T^^K~'^^w± (A), where 
w± (A) are entire i/i-periodic functions. However, using the asymptotic behavior of q and 
we get that w± (A) are bounded at infinity in the strip 10" (A) | < h/2, and hence on C. They are 
thus constant. This proves the uniqueness of solutions for a given choice of r^'s and hence 5fc's. 
Indeed, up to a normalization constant, any solution q satisfying to the requirements stated in 
point (ii), is of the form 



Ntt 



g+ (A) - eg,- (A) 



g(A)=e^" ';^ '"^ (C.U) 
n sinh T (-^ - ^k) 

k=l " 



As the solution q is entire, it has a vanishing residue at A = 5k, k = 1 . . . , N. Therefore, 
the quantization conditions for the Toda chain appear as the set of — 1 conditions that 
q (A) has a vanishing residue at 5k, k = 1, . . . , N supplemented with the N^^ quantization 
condition for the overall momentum : J2p=i '^p — Yl,p=i — ^- Note that it follows from 
^ [Qt^ Qj] i^k + inh) = that if q has a vanishing residue at a 5k then it also has a vanishing 
residue at 5k + ink, n E Z. Therefore, there is indeed only a finite number of constraints of 
the parameters 5. The explicit form of these quantization conditions is then indeed as given in 
(EH). 

Conversely, if the quantization conditions are satisfied, then by taking q as in (IC.llI) . one obtains 
an entire solution with the desired asymptotics. □ 
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